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Abstract. Using a method of stochastic perturbation of a Langevin system 
associated with the non-viscous Burgers equation we construct a solution to the 
Riemann problem for the pressureless gas dynamics describing sticky particles. 
As a bridging step we consider a medium consisting of noninteracting particles. 
We analyze the difference in the behavior of discontinuous solutions for these 
two models and the relations between them. In our framework we obtain 
a unique entropy solution to the Riemann problem in ID case. Moreover, 
we describe how starting from smooth data a 6 - singularity arises in one 
component of the solution. 

Introduction 

We propose a method for solving the Riemann problem as well as for describing 
the formation of singularities for the pressureless gas dynamics system and a natural 
extension of it. The system of pressureless gas dynamics is very important since 
it is believed to be the simplest model describing the formation of structures in 
the universe (e.g. [27]) and plays a significant role in the theory of cooling gases and 
granular materials [7] . It is a system consisting of two equations for the components 
of the density / and velocity u expressing the conservation of mass and momentum 

d t f + div x (fu) = 0, (1) 

$(/«)+ V B (/u®u) = 0. (2) 
It first appears to be very simple, however a closer analysis reveals that it has some 
peculiar features due to its non strict hyperbolicity. The system has attracted a sig- 
nificant interest in the last decades and has been investigated quite intensively. In 
particular, it is well known that the arising in the velocity component of unbounded 
space derivatives implies the generation of a 5 - singularity in the component of the 
density. Therefore for this system one needs to define a generalized or measure- 
valued solution of a special kind. This was done in [2D], 0, [113 , [35] , [2j|] , [23], [S], 
where the authors used different techniques (vanishing viscosity, weak asymptotics, 
variational principle, duality) to define the solution and prove its existence. Fur- 
ther, the Riemann problem for the pressureless gas dynamics was studied (e.g. [29] . 
[16|.[28|). including a singular Riemann problem with a <5 - singularity concentrated 
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at the jump at the initial moment. Nevertheless, even in the ID case there are 
certain open problems, not to mention those present in the higher dimensional sit- 
uation. In particular, there is a problem concerning the uniqueness of solutions. 
Both in the case of rarefaction and contraction it is possible to construct a whole 
family of solutions to the Riemann problem satisfying the integral identities and 
entropy conditions that are used to single out the unique solution in the strictly hy- 
perbolic case (see [12], [19] for details). Further, the process of singularity formation 
was described up to now only in a very special situation f|13j.[Tl]). 

Our method allows to find a unique solution to the Riemann problem with ar- 
bitrary smooth left and right states. The most clear and explicit results we get 
concern the case of constant left and right states. We restrict ourselves to the ID 
case, however the formulas that we use are written in a similar way in any dimen- 
sions, and our technique can be straightforwardly extended to the case of higher 
dimensions. Of course the situation in higher dimension is much more complicated, 
nevertheless, there are no fundamental obstacles to be faced when applying our 
method. 

Further, our method is constructive. It allows to describe the behavior of the 
system starting from any smooth initial data. In particular, it is possible to describe 
the singularity formation including the time, position and value of the amplitude 
of the 8 - function in the component of the density basing on the initial data. 
Moreover, it is possible to describe the behavior of a solution after the critical time, 
in particular, the position and the amplitude of the 8 - singularity. 

Let us describe shortly the method in the ID case. First of all it is evident that till 
the solution to the system of pressureless gas dynamics keeps its C 1 - smoothness, 
the velocity satisfies the non-viscous Burgers (or free transport) equation. We 
write it in the Langevin form and introduce a stochastic perturbation along the 
trajectories of the particles. Further, we consider the position x and velocity u of 
the particles as random variables and find the common probability density P(t, x, u) 
in the space of positions and velocities as a solution to the corresponding Fokker- 
Planck equation. Then we introduce a pair of variables: the density of particles 
p(t, x) and the conditional expectation of velocities u(t, x) at fixed coordinates (and 
time) . 

We define a generalized solution in the sense of free particles (FP), a pair 
((fpp(t, x),UFp(t, x)), as a special double limit of (p,u) as the parameter a of 
stochastic perturbation and the parameter s of the approximation of initial data 
go to zero. 

We prove that the pair (p, u) satisfies a gas dynamics system with viscous and 
integral terms. The viscous term is the usual one for the viscous approximation of 
the solution and it vanishes as the parameter of stochastic perturbation a goes to 
zero. The integral term vanishes if upp, a part of the FP - generalized solution, 
is continuous and persists otherwise. Thus, the FP - generalized solution to the 
pressureless gas dynamics system solves in fact the extended gas dynamics system 
with an integral term in the usual sense of integral identities. This integral term 
can be considered as a spurious pressure, it is equal to the dispersion of u with 
respect to P x (t,x,u) (P x denoting the derivation with respect to x). 

The FP - solution corresponds to the model describing the behavior of a medium 
consisting in a micro-level of non-interacting particles. In the component of density 
of the FP-solution the 8 - singularity can arise only from the domain where of the 
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initial velocity has a special form (see Sec. 6). The S - singularity does not arise 
in the FP- solution to the Ricmann problem with constant left and right states 
provided the S - singularity was not concentrated at the jump initially: instead of 
the S - singularity we have an overlapping domain of a non-zero measure (a spurious 
pressure arises in this overlapping domain) . The FP-solution is interesting in itself, 
moreover, it can be used to construct a solution to the sticky particles models, 
where the particles are assumed to move together when they meet (the name of 
sticky particles model is used as a synonym of the pressureless gas dynamics). We 
propose a method based on the conservation of the mass and momentum that allows 
to reduce the FP-solution to the solution of the sticky particles model and find the 
position and the amplitude of the S - singularity in the density component. The 
spurious pressure degenerates in the process of above reduction. It fact, the problem 
of collapsing the overlapping domain into a point only arises for the compression 
waves, since for the rarefaction waves the initial jump in the velocity decays into a 
smooth profile, such that the FP-solution and the solution in the sense of integral 
identities coincide. 

The paper is organized as follows. In Sec.l we consider the model of motion for 
free particles perturbed along their trajectories and introduce the integral charac- 
teristics of the medium consisting of these particles, in particular, the mean velocity 
at a fixed coordinate u. In Sec. 2 we study the properties of u and prove that the 
limit of this value as the parameter of the stochastic perturbation go to zero (pro- 
vided it is smooth) takes part of solution to the pressureless gas dynamics system. 
In Sec. 3 we give a notion of generalized solution in the sense of free particles (FP) 
to the Cauchy problem for the pressureless gas dynamics system. In Sees. 4 and 5 in 
ID case we find the FP-solution for the classical and singular Riemann problems, 
respectively. In Sec. 6 we describe the arising of singularity from smooth initial 
data for the pressureless gas dynamics system. In Sec. 7 we discuss the difference 
between the FP-solution and the generalized solution in the sense of integral iden- 
tities. In Sec. 8 we propose a method of changing the FP-solution to the solution 
in the sense of integral identities, and discuss the solution to the Riemann problem 
with non-constant states and the evolution of the singularity arising from smooth 
data. In Sec. 9 we extend the method to a class of systems that can be obtained 
from the scalar conservation law with a convex flux. Sec. 10 is a conclusion where 
we discuss the related approaches and methods. 

1. Stochastic perturbation of the Burgers equation 

Let us consider the Cauchy problem for the non- viscous Burgers equation: 

d t u + (t»,V)« = 0, t > 0, u(x,0) = u (x) GC 1 (M™)nC 6 (M n ). (3) 

Here u(x,t) = {u±, ...,u n )(x,t) is a vector-function IR n+1 — > R™. 

It is well known that solving this equation in the smooth setting is equivalent to 
solving the system of ODEs 



for the characteristics x = x(t). 

We associate with ((4]) the following system of stochastic differential equations: 



±{t) 



u(t,x(t)),t>0, ii(t,x(t))=0 



(4) 



dX k (t) = U k {t)dt + crd{W k ) 
dU k (t)=0, k = l,...,n, 



(5) 
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X(0) = X, U(0) = U, 

where X(t) and U(t) are considered as random variables with given initial distri- 
butions, (X(t), U(t)) runs in the phase space ffi™ x M. n , a is a real strictly positive 
constant and (W)t — (Wk)t, k — 1, ...,n is an n - dimensional Brownian motion, 
t> 0. 

System (3) describes the free motion of a particle that "does not feel" the other 
particles. We assume that initially at time t = these non- interacting particles are 
distributed with a density f(x) and denote by P(t,x,u) the probability density in 
position and velocity space for the solutions of {5}, x € K™, u € W 1 at time i. 

The stochastic system described by ([5]) is associated with the deterministic non- 
viscous Burgers equation ([3} in the sense that the deterministic characteristics 
(j4| are replaced by a stochastic perturbation of the characteristics as described 
by ([5]). System ([5} with this interpretation is what we understand as "stochastic 
perturbation of the Burgers equation." 

Let us introduce the function 

J uP(t, x, u)du 

^ X) = m }p(t,x,u )d u > ^ ^ (6) 

This value @ can be interpreted as the conditional expectation of U for fixed 
position X [TD] . If we choose as initial distribution 



Po(x,u) = S(u - u (x))f (x) = Y[ 5{u k - (u (x)) k )f Q (x), 



(7) 



k=l 



where fo is an arbitrary sufficiently regular nonnegative function such that J fo(x)dx 

R" 

1, then u(0,x) = uq(x). Certain properties of u(t,x) have been established in pQ( 
see also [2] for another type of stochastic perturbation) . 

The density P — P(t, x, u) obeys the Fokker-Planck equation 



dP 



u '- 



d 



2 a dxl 



(8) 



subject to the initial data ([7]). 

We apply the Fourier transform to P(t,x,u) in ©, 1(7} with respect to the 
variables x and u and obtain the Cauchy problem for the Fourier transform P = 
P(t,\,£) of P(t,x,u): 



P(0, A, = / e- 4 ^e- 4 «^» jo( S )ds, A, f € 



(9) 
(10) 



Equation ((5} can easily be integrated and we obtain the solution given by the 
following formula: 

P(t, A, = P(0, A, £ + \t)e-i° 2 M 2t . (11) 
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The inverse Fourier transform (in the distributional sense) allows to find the 
density P(t,x,u), t > 0: 

R™ R" 

1 ' / " e i(A,*) e i«,«) f e -i(A, fl ) e -i«+At, U0 («)) f Q ( s )dse-^ 2 ^ 2t d\d( = 



(2tt 



2)i 



(27T) 



■ •• ' O i(?,u-«o(s))jt / -?<r f ( A ^ J " 



Ma) / e l ^"-"°WJ<i£ / V ^ V ^ dXds 



' ^ <5(u- uo(s))/o(s)e _ ' l ' 0<S 2^; x ' ds, i > 0, i £ 1". (12) 



(Verier)™ 

Then we substitute P(t,x,u) in ^ and get the following expression for u(t,x) 
(sometimes we insert a label a to stress the dependence on this parameter): 

ti (s)t + s-a;| 2 

/ u (s)f a (s)e ^ ds 

u(t,x) =u a {t,x) = — i„ 0(a )t+.-x|=! • ( 13 ) 

/ fo(s)e ds 

Remark 1. The integrals in I113\) are defined also for a wider class of fo than the 
probability density of the particle positions in the space at the initial moment of 
time. If the integral J fo(x)dx diverges (for example, for fo(x) — const), we can 

consider the domain [—L,L] n , where L > and use another definition of u a {t,x): 

J u (s)f (s)e 2,^* ds 

U (T (t,x) = lim — |„ n(s)t+s -tt|2 ( 14 ) 

J Jo(s)e ^ ds 

[-L,L] n 

(provided the limit exists). Evidently, this definition coincides with H3\) for fo G 
Li(K"). 

2. Properties of velocity averaged at a fixed coordinate 

The following property of u(t, x) holds: 

Proposition 1. Let uo and fo > be functions of class C 1 (R") n C/,(IR™). If 
t*(uo) > is a moment of time such that the solution to the Cauchy problem f3|) 
with the initial condition uq keeps this smoothness for < t < t*(uo) < +oo, then 
Ucrit, x) tends to a solution of problem (0) as a —>■ for any fixed (t, x) G R" +1 , < 
t < t*(uo). 

Proof. Let us denote by J(uo(x)) the Jacobian matrix of the map x i — ► uq(x). 
As it was shown in [25 (Theorem 1), if J(uo(x)) has at least one eigenvalue which 
is negative for a certain point x G W 1 , then the classical solution to ([3|) fails to 
exist beyond a positive time f*(uo)- Otherwise, t*(uo) — oo. The matrix C(t,x) = 
(I + tJ(uo(x) j), where I is the identity matrix, fails to be invertible for t — t*(uo). 
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The formula (I13| (or l|14[)) implies, using the weak convergence of measures and 
the fact that fo and uo are continuous and bounded 

/ uo00/o(*) Hm (v ^ CT) „ e-'"° ( tt; 3=1 rfs 

,. « ,. x R" CT u v 1 

J uo(s)fo(s)S p ^ tX>s )ds 

with s) = uo(s)t + s — x , where <5j, is the Dirac measure as y € M.". We can 

then on the basis on the invertibility of C(t,x) use locally the implicit function 
theorem and find s = St, x (p). Therefore, 

/ Uo(st,x(p))fo(st,x(p)) det(C(t, st^ip)))^ 1 5 P (ds t ^ x ) 

lim u a (t,x) = — ; r = u (s t ,x(0))- 

CT ^o J fo(st,x(p)) det{C(t,st,x(p))) 1 Sp{dst, x ) 

R" 

Let us introduce the new notation so(t,x) = St )X (0). Then the following vectorial 
equation holds: 

uo(so(t, x))t + so(t, x) — x = 0. (15) 
Let us show that u{t,x) — uo(so(t, x)) satisfies the Burgers equation, that is 

n n 

dj(uo t i)(s j) t + uojdk(uo t i)(s 0t k)x j = 0, i = l,...,n, (16) 
j=i 3 ,k=i 

and uo(so(0, a;)) = Uq(x). Here we denote by ito,i and So,j the z - th components of 
vectors Mo and sq, respectively. 

We differentiate (JT5]) with respect to t and Xj to get the matrix equations: 

n 

C»j (go,j)t + uq,% = 0; i = l, ...,n, 

i=i 

and 

n 

Cifc (so.fc)^ + % = 0, i,j = 1, n, 

fe=i 

where is the Kronecker symbol. The equations imply 

n 

(so,j)t = - ^2 (C _1 )ij M0,i, (S0,k)x 3 = -(C -1 )jfc. (17) 
i=l 

It remains now only to substitute ([TT)) into (fTfJ|) . 
Further, (fTS"]) implies uo(so(0, a;)) = tto(x). □ 

It is important to note that So(t,x) is unique for all t for which the solution to 
the Burgers equation u(t, x) is smooth. 



Remark 2. Proposition^ can naturally be extended to the class of functions fo 
such that there exists a sequence /g £ C 1 (M™) n C&(K n ) converging to fo as e — ► 
almost everywhere. In this case u a (t,x) tends to a solution of problem |3|) as er — ► 
almost everywhere on (t,x) € R n+1 , < t < £*(wo)- 
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Let us set p(t, x) = J P(t, x, u)du, t>0 7 xE M™. From (fl~2|) we have 

1 



p(t,x) = p a (t,x) 



r 2irtay 



\u Mi+»-»r 
fo(s)e ^ rfs, t > 0, x S 



(18) 



Proposition 2. The scalar function p(t, x) and the vector- function u(t, x) = (ui,..,u n ) 
defined in (0|) solve the following system: 



- + d iv x (pu) = -a ^_ 



(19) 



fe=i 



d(puj) 
dt 



\7{pu lU ) = -(T 2 ^ 



d 2 (pui) 



2~ ^ dz 2 . 



(ui - Mi) ((u - u),V x P(t,x,u)) du, 

(20) 



i = 1, ■•, n, t > 0. 



Proof. The equation (fl9| follows from the Fokker-Planck equation (J8J) directly. To 
prove (|20p we note that the definitions of u(t, x) and p(i, x) imply 

— — I uP(t,x,u)du= / uP t (t, x,u) du — 



dt dt 



U(U, V x P(t, X, u))du + -cr 2 ^2 -q^I' 

k—1 ^ 



(21) 



where P t = J^-P. 



Further, for i = 1, .., n we have 



Mfc P du 



' J Uk P du'' 



Uk P du 



dx k 



\ 



J Pdu 



= / in uk P Xk du 



J Ui P Xk du J P du — J Ui P du J P Xk du 
m. i',!,u- — 7 ^ = (22) 



J Pdu 

v J 

{u k Ui + u l u k - u k Ui) P Xk du, i,k = l, n, 



_ d 



P. 



with P r , - - - — 

• Lk ox k 

Equation (J2(5]) follows immediately from (21]) and [[22]), □ 
In the one-dimensional case (fH?|) - (l2TH) has the form 

d t p + d x {pu) = ^<r 2 d 2 x p, 

d t {pu) + d x (pu 2 ) = ^cr 2 d xx (pu) - / (u - u) 2 P x (t,x,u)du. 



(23) 
(24) 



Let us set x) = lim x) and x) — lim u{t, x). 



8 



ALBEVERIO, KORSHUNOVA, ROZANOVA 



Proposition 3. Assume that (f(t, x), u(t, x)), the limits of (p, u) as a — > 0, are C 1 
- smooth bounded functions for (t,x) € O := [0,i*(wo)) x R™, t*(uo) < oo. Then 
they solve in Q the pressureless gas dynamics system (QP, @). 

Proof. As follows from Proposition [1] the function u(t,x) is a C 1 - solution of the 
non-viscous Burgers equation. Further, (|23|) is a linear parabolic equation with 
respect to p, hence the limit as a — > reduces it to the continuity equation |T]). 
Equation ([2|) is a corollary of the non- viscous Burgers equation and (pj. □ 

Remark 3. Proposition^ implies that the integral term on the right-hand side of 
I120\) vanishes as a — > in the case of smooth limit functions f and u. Let us prove 
this fact alternatively. Indeed, as follows from we have as a — > 



(V27rfcr)' 

JK" 

1 f | UQ ( 3 )t +a - 



(ui — Ui) ((it — u),V x P(t, x,u)j du 



_ \UQ (S)t+S-X\ 2 

fo(s) (u 0i (s) - Ui) ((uo(s) - u) , \7 x e ) ds 



(\T2irta) n to- 2 



fo(s) (uoi(s) - Ui) ((u (s) - u) , uo(s)t+s-x) ds = 



' e / (s) (u 0l (s) - Ui) ((tto(s) - u) , u (s)t+s—x) ds +o(a), 



{^Thrio-yto- 2 

where U So u,x)( e ) is an e - neighborhood of the point so(t, x) (see the proof of Propo- 
sition^), £=1, ..,n. Further, 

\u 0i (s) - Ui{t,x)\ = |(ttof(s) - Ui(t,x)) + (ui(t,x) - Ui(t,x))\ = 

|(ttoi(s) - u 0i (s (t,x))) + (u 0i (s Q (t,x)) - Ui(t,x))\. 

For every fixed x andt € (0, t*) and for every a > there exists s(a) > such that if 
s 6 U So n x \(s), then \(uoi(s) — uoi(so(t, x))) < a. Moreover, since u ai (t, x) — > w(i, x) 
as o~x — > (we Ziaue renamed the parameter), then for every a > f/iere exists 
o~i(a) > swc/i that for sufficiently small a \ we have \u,i(t,x) — Uoi(so(t,x)\ < a. 
Thus, 

' e 2^ t / (s) (tt Qi (s) - Hi) ((tto(s) - u) , uo(s)i+s-x) ds| < 



(V27rtcr)™to 2 

C/ S0(t ,x)(e) 

Const • — — — f e /o(s) \ u o( s )t + s — ^1 ds, 



l^nta) r 

where the constant does not depend of a. The latter integral tends to zero as a — > 0. 

In fact, to prove that the integral term vanishes as a — » 0, we have used only the 
continuity of u and the boundedness of fa. 

However, as we will show in Sec^ if we put instead of u a discontinuous func- 
tion, this integral term does not vanish. 
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3. Generalized solution in the sense of free particles 

Being inspired by the fact that the formula (fl3|) makes sense also for discontin- 
uous initial data (/o(x), uq(x)), we give the following definition for any dimension. 

Definition 1. We call the couple of functions (fpp{t, x),UFp(t, x)) a generalized 
solution to the Cauchy problem in the sense of free particles (FP- generalized 

solution) subject to initial data (/o(x), Mo(x)) £ L^ oc (M") nLoo(R™), if for almost 
all (t, x) £ M+ x R™ 

f F p(t,x) = lim (lim p e a (t,x)), 

e—>0 <t— »0 

UFp(t,x) = lim(lim u e a (t,x)), 
where (p £ a {t,x),u e a {t,x)) satisfy the system \2S\) , \2I$ with initial data 

fo = Ve * fo, "o = Ve * Uq 

where n e (x) is the standard averaging kernel. 

Remark 4. The properties of the standard averaging kernel |21j imply that /q and 
Uq) belong to the class C°° and 

Km fo (x) = fo (x) , lim vf (x) = u (x) , 
for almost all fixed x £ R" . 

Remark 5. As we will see below, if ufp is continuous, the FP- generalized solution 
is a solution to 0), for example, in the sense of integral identities. However, 
if ufp is discontinuous, the FP-solution solves a different system, namely one that 
differs from ^j), @) by an integral term corresponding to a spurious pressure. Nev- 
ertheless, using the FP-solution we can solve (QP, (0j itself. 

Definition 2. We call the pair (/o,Mo) a monotonic approximation of initial data 
(fo(x),u (x)) e ]Lf oc (R™) n L 00 (R n ), if 

• /o an d u o are from the class Cf,(R™), moreover, they are from C (R n ) 
almost everywhere; 

lim fo (x) = fr (x) , lim u E (x) = u (x) , 

for almost all fixed x £ R ra ; 

• for sufficiently small £ and almost all fixed (x, t) £ R™ +1 every root Sk of 
the equation Uq(s) t + s — x = belongs to the neighborhood U Sk (e) of the 
root Sk of the equation Uq(s) t + s — x = 0. 

Proposition 4. The FP - solution (fFp(t,x),UFp(t,x)) to the Cauchy problem 
(0) with initial data (fo(x),u (x)) £ L| oc (R n ) nL. 0O (M n ) does not depend of 
the choice of the monotonic approximation (/q,Mq). 

Proof Let us choose two the monotonic approximations (/qi, u§i) and (/q 2 , Wq 2 ) 
such that 

!im / £ i(x) = lim / e 2 (x) = / (x), 

£—►0 £— >0 

lim Uoi(x) = lim Uq 2 (x) = Uo(x) 

e— >0 e— >0 

for almost all fixed x £ R". Then the couple 

(foi u o) = (foi ~ /o2j u 01 ~ u 02) 
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can be considered as initial data for the problem ([I])-®. To prove the proposition 
we have to show that the respective solution is identically zero almost everywhere. 
Indeed, from (JT5J) we have for any t > 0, and almost all x G M. n 

f FP (t, x)=J2 ( / fo(s)S(s - sU(t, x))ds 
Hm ( / r 02 {s) (8(s - *?,*(*,*)) - S(s - sl^x))) ds 



Y,^foi(4A t ^))~fo2(4.k(t,x))) - lim(/5 ! (4 fc (*,x))-/^( a ? |fc (t, a; )))]=0. 

' e— >0 e— »0 

fe 

Here s| fc (t, x) is the fe-th solution (fe = 1, 2, if, i = 1, 2) of equation 

ul i (s)t + s-x = 0, i = l,2. 

We have used the fact that |sf ie) — s| x)| — > 0, fe = 1,2, if as e — > 0. 
Analogously proceeding from (113|) . we prove that upp(t, x) = for almost all 
t > 0, a: G R™. □ 



4. The classical Riemann problem in the FP sense for the ID case 

For the sake of simplicity we restrict ourselves to the one-dimensional case and 
consider the following initial data: 

f (x)=f 1 (x)+6(x)f 2 (x), (25) 

u (x) = ui(x) + 9(x)u 2 (x), (26) 

where 9 is the Heaviside function with jump at xo (without loss of generality we 
shall assume xq = 0), u\, u 2l fi, f% are continuously differentiable functions. We 
shall dwell first on the case where /i, f2,Ui, u 2 are real constants. 

According to Definition [T] wc must consider the smoothed initial data instead 
of (j25|) and (f26|) . It follows from Proposition [1] that we can choose any couple 
of smoothed initial data. It will be convenient to consider the piecewise linear 
monotonic approximation of initial data of the form 

fi, x < -e, 

/u<<> = { Y s x + fl + f' - £<x<£ > W 

Ji + fi, x>e, 
ui, x < — e, 

Uo(x) = { 7^; x + u i + y > -^<x<e, (28) 

Ui +U 2 , X > £, 

where fx, f 2 , u± and u 2 are real constants. 

From psp we can find the density p e a (t, x) corresponding to the smoothed initial 
data (/q (x), Uq(x)) (below we shall omit the index a for short). 

Let us set: 



Ei(t, x, s) = exp 



(s — x + U\t) 
2^H 



(29) 
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E 2 (t, X, s) = exp 
E%{t, x, s) = exp 



(s - x + (ui + u 2 )t) 2 



2aH 



2^H 



(30) 
(31) 



Then 



P £ (t,x) = 



1 



/i£i(f, x, s)ds + / {fx + f 2 )E 2 (t, x, s)ds+ 



I (^s + f 1 + ^)E 3 (t,x,s)ds 



We denote the third integral by and set p = 
We then have: 



C^ + l)s+(u 1 + ^)t-x 



I E > a = L s (t,x)a / pe~-dp + F e {t,x)a / e^-dp 



L e (t, x)a ( e^^^~ - e~^ r ) + F £ (t, x 
1 



C'_ /cry/t 
CI 



' + 
(T 



- 3> 



6*1 



(32) 



where $(a) 



2tt. 

(mi + u 2 )t — x + s, and 



F e (t,x) 



J e 2 dx is the Gauss function, Cf. = itit — a; — e, 



L e (i,a;) = - 



2i/ 2 e 



V^(u 2 i + 2e) 2 ' 



2s 



u 2 t + 2s 



(33) 



It can easily be seen that lim F £ (t, x) = 0. 

e— >0 

Finally, we get 

f(t, x) = A$ (^] + (A + / 3 )* f — % ) + I{' 



ry/ij 



(34) 



To find u(t, x) we compute the numerator in formula |T 



1 



y/2nta 
where 



«o( s )/o( s ) e ~ ° ^ ds = u lP £ (t, x) + u 2 {fi + / 2 )$ ( + L 



Ci 



2 i 



je,a = G e >°(t,x)+K e {t,x)<j\e 



(01)-= (ci)- 



2tcr- g 2tCT 



(35) 
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N £ (t,x) 
Here we used the notation 



K s {t,x) = - ! =^^— T F e (t 1 x) + J^Lu £ (t,x), 



2Tr(u 2 t + 2e) 



2V2ne 



N £ (t,x) 



We recall that F £ (t, x) was determined by (|33|) . 
It is easy to deduce that lim N £ (t, x) = 0. 

Thus, we have the following result: 



u £ {t, x) = Ux + 



u 2 ( 



(36) 



CI 

where 7^' cr and /| :CT are given by (f3"2|) and ([33]) , respectively. Note that — > 

(7 

— > ±oo as cr — > 0. 

It can be checked that the initial conditions are satisfied, namely p(0, x) = fo(x) 
and u e (0, x) — u*q(x). 

Now we can find the generalized solution to the Riemann problem in the form: 

f FP (t,x) = lim (lim p £ (t,x)), 

e^O cr— »0 



upp(t,x) = lim( lim u e (i, x)). 

Let us introduce the points if = u\t — e and x| = (lii +U2)t + £. Their velocities 
are u\ and u% + u%, respectively. 



We consider two cases: 

1. U2 > (velocity of the point x\ is higher than the velocity of the point xf) 
At first, we can find f £ (t, x) = lim p £ (t, x) from (f3"l)) . It is easy to see that 

(7 — >0 



A, 

F s (i } x), 

2 — +2 F ( t > x >> x 

h + /2, 



Let us note that this formula contains F £ (t,x) and \im F £ (t, x) = 0. Thus, we 
obtain the following result for fpp(t,x) = lim f e (t,x): 

e— +0 
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A, 
l± 

2 ' 

f F p(t,x) = { 0, 



X < Uit, 
X = Uit, 

U\t < X < (Ul + U2)t, 

x = (ui + u 2 )t, 
fl +/2, x > (ux + u 2 )t. 



h + h 



Further, from (|36[) we find the solution of the gas dynamic system with smooth 
initial data u e (t,x) = lim u e (t,x) as follows: 

er— >0 



u £ (t, x) 



It can be shown that 

,. N e {t,x) 
hm ; 

F e {t,x) 



lim 



Ul, 



"I 



N £ (t,x) 
F £ (t,x)' 

Ui + U2, 



x <C x^ , 

3/ < ~^~ X 2 ^ 



"2 



(a; - (in + — )t) + — 
ou 2 t + 2e y y 1 2 ' ' 2 



Thus, we get the following solution for u(t, x) = lim u £ (t, x) 

e— >o 



Ui. 



Ui, X < Uit, 

u FP (t, x) — { -, Uit < x < (ui + u 2 )t, 

Ul + u 2 , x > (ui + u 2 )t. 

We can see that the velocity includes the rarefaction wave (see FigfT]). This is a 
well known self-similar solution to the Riemann problem with constant left-hand 
and right-hand states for the Burgers equation (|29|). 




U t (UjtUjIt 



Figure 1. Density and velocity, u 2 > 0. 
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It is interesting to note that if we first compute the limit in e we get the solu- 
tion u(t, x) = ui + u 2 8(x — (u\ H — — )t), which is unstable with respect to small 
perturbations. 

2. u 2 < (the velocity of x\ is higher than the velocity of if ). From (I34|) and 
((36)1 we find as before (see FigfJ]) 



f FP (t,x) = < 



u FP (t,x) 



fl, 

3/i + h 
2 ' 
2/1+/2, 
3/i + 2/ 2 
2 

/1 + /2, 




/1 + /2 
2/i + h 
Ui + u 2 , 



X < (Ui + U2)t, 

x = (ui + u 2 )t, 

(ui + u 2 )t < X < U\t, 

X = Ult, 
X > Ult, 

X < (ui+ u 2 )t, 

u 2 , (ui + U 2 )t < X < Ult, 
X > Ult. 



2 f x +f 2 

♦ 

f, +fj 



lu 1 + u z ) X, 



Figure 2. Density and velocity, u 2 < 0. 



5. The singular Riemann problem in the FP sense 

By the singular Riemann problem we mean the Cauchy problem with the fol- 
lowing data: 

fo (x) = fi + f 2 0(x - x ) + f 3 S(x -x ) = f reg + f sing , (37) 
uq(x) = ui + u 2 9(x — xo), (38) 
which differs from the classical Riemann problem (pj)) . by a singular part fg ing , 
the S - function of constant amplitude f 3 concentrated at the jump of the density. 
We again set xo equal to zero (without loss of generality). 

Thus, we should adapt the definition of the generalized solution to the case 
of a singular density. As before we want to approximate initial data by smooth 
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functions, the 8- function will be instead naturally approximated in the space 2?'(R) 
of distributions. 

Definition 3. We call the couple of functions {fFp{t,x),UFp{t,x)) a generalized 
solution to the problem (QP , ifl?[J, {Sty , i38\) in the sense of free particles (FP), if 

f F p{t,x) = hWlim i p% reg (t,x)) + lim(limp £ i (t,x)) = f reg (t,x) + f sing (t,x), 

upp(t,x) = lim(lim u s T (t,x)), 

where the limits are meant as pointwise at almost all points except for the singular 
part 

lim(lim/£ H ng (t,x)), 

where the outer limit in e is meant in the sense of distributions. 

Here (p% ireg {t, x) + p £ a sing (t, x), it e (t, x)) satisfy the system l23\) . \2J$ with initial 
data (fo reg (x) + /g sing {x), Uq(x)) from the class C 1 (M n ) such that 

: e— >0 



for almost all fixed x £ 1" and 



the limit in e being in the sense of distributions. 



We are going to solve the singular Riemann problem with constant left and right 
states. To approximate the 5-function we use the well known fact that 
1 



e 2e 



S(x), 



inX>'( 



The part of the solution that relates to the regular part of the density (/3 = 0) is 
found in SecHJ Now we have to calculate the singular part of the density p% sing (t, x) 
1 

for /q sing {x) — fa e 2 e . As before we omit the index a. 



%gK ' - \/2lTE 

From II18D we obtain: 



+ exp 

h 

^2ir(c-H + e) 



exp 



y/2Tr(an + e) 
((tti + u 2 )t - x) 2 



(uit — x) z 



2(<7 2 t + e) 



2{a 2 t + e) 
D 



Gl' e {t,x)<5> 



Dl> 



where 



(uit - x)^/e 
^Ja 2 t + e 



aVt 
Gr{t,x)>S> 



- y/e(a 2 t + e), 



Do 



rVi 



u 2 — 



((mi +u 2 )t~x)^/e i 2 ,, , 

+ \j£(p H + e), 

Vcr^t + e 



and Jn 



I ex P 



s 

2e 



(s~x + (^s + Ul + ^)t) 2 
2^t 



ds . 



+ 



(39) 
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It can be checked that 

je.cr 



2 y/2^VA^ 

where 



exp 



2 i 



2A a ' E 



- $ 



K°_ 



o\ft 



B° 



le 

:= £ (^ + l)((u 1 + ^)t-x) 



K ± e = "4= (VA°*l±e) + B a - 

\ E \ 



(40) 
(41) 
(42) 



It is easy to see that 

D\ = lim D^' e = u\t — x — e, D\ = lim D%' e — {u\ + u 2 )t — x — e, 

Kl = = + l)(±e) + («i + f )t - x, 

i.e. iff. = £)f = Cf , iff. = Df = CI. Therefore, in this case we also have two 
jump points if = U]f — e and i| — ( u i + u 2)£ + £• 

We have to consider two cases it2 > and U2 < 0, as before. But for fo in9 (x) = 
= fz <5(ic) we obtain the same result: 

fsing{t,x) = lim lim p E si ng (t,x) = ^f 3 (S{x - u x t) + S(x - (m +u 2 )t)). 



1 



Thus, we have two 5-functions with equal amplitudes —F(t), which move with the 
jump points. 

Further, we can find the regular part of the density f reQ (t, x) — lim lim p e rea (t, x) 
where 

pt(t, x) = /i$ (^-] + (h + / 2 )$ (-^) + h- 



(43) 



Analogously, from (fT3"|) we can calculate the velocity u(t,x): 



"2 



u E (t, x) = u%+- 
where 



(/l + / 2 )$ (- 



C: 



_+_ ) J- 7^°" 



p £ (t,z) 



jsi<t y*3 

V Z7T 



1 /«a w 2 S° 



V 2 2eA a 
G a ' £ (t,x) = exp 

Gf' e (i,x) = exp 



U2V te 
2eA^' 



e 2t ° 



if! 



- $ 



((ui + ^t-s) 



e 2t " 



kit 

<j\fi 

2 1 



((ui + u 2 )t — x) s 
2(aH + e) 
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The expressions for p £ (t,x) — p £ eg (t,x) + p E a i ng (t,x) and I 2 a have been given 
in (|43l), (091) and (35]), respectively. 
1. i*2 > then: 



u £ (t,x) = < 



Ui 



T £ G £ (t,x)+N £ 



fa 



2s 



/2~ne 2e + u 2 t 



G £ (t,x)+F £ 



where 



T £ = lim 



u\ + u 2 , 



1 (u 2 u 2 B° 



2eA i7 ' s 



Therefore 



UFp{t,x) = \\mu £ {t, x) = < 

wi + u 2l x > X^Ul + u 2 )t 
Further, from P3")) , (|3T)|) we find the density as follows: 



Ul, 

a; 



x <C ai^ , 

^1 ^ ^ *^2 ' 
X > X 2 , 

u 2 2e((ui + y)t-a;)M2 
Y (2e + u 2 i) 2 Ve 



a: < uit, 

< a: < (iti + u 2 )t, 



f F p{t } x) = -f 3 (5(x-u 1 t) + 5{x-(u 1 +u 2 )t))+ < 



2. If u 2 < 0, we have 



fi- 
ll 

2 ' 
0, 

h + h 



X < Uit, 
X = Uit, 

Uit < X < (ill + u 2 )t, 

x = (iti + u 2 )i, 



/l +/2, a; > («i + u 2 )t, 



u e (t, x) = < 



Ul + 



x < x\ 



Mh + h) -N £ + ^=(G| - sT £ G £ ) 

VZTTE 



where G\ = exp 
Thus, 



Vx + f2 + F £ 
U\ + u 2 , 

(uit — x) 2 



2f ' 
--(G\ + G £ 2 -——-G £ ) 



I2me 



2e + u 2 t 



X > x 2 , 



2e 



, G| = exp 



((ui + u 2 )t - x) 



2 1 



2e 



Ml, 



X < (ill + M2)^, 



u FP (t,x) — ^ ui + ^ ~^/^) U 2 , (Ul + U 2 )t < X < Uit, 



2/i + / 2 

,1*1 + M 2 , 



a; > uit. 
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As a consequence, for the velocity we obtain the same result as in the non-singular 
case. Analogously, 



fx, x < {ui + u 2 )t, 

3/i + h . ^ s, 
- , x={u l +u 2 )t, 

2/i + f 2 , (ui + u 2 )t < x < u\t, 



X = U\t 



3/i + 2/ a 
2 

/1 + /2, x > U\t. 

6. Singularity arising from smooth data 

We are going to show that at the point of formation of a singularity from smooth 
initial data in the solution to the pressureless gas dynamics model a S - function 
appears in the density component. For the sake of simplicity we again restrict 
ourselves to the ID case. 

Theorem 1. (Asymptotics of the approximating solution) Let the initial data 
(/o; uo) for the system (QJ), @) be at least C m - smooth and bounded, m > 2. 

Assume that there exists an instant < t* < 00, such that t* = inf I 77— r ) 

and u( k \s*) = 0, k = l,...,m— 1, however it( m )(s*) does not vanish at the point 
s*(t*, x*), where s* is a solution to the equation Uq(s) t* = x* — s. Here x* is such 
that the line y = x " t ~ s intersects the graph of the initial velocity y = uq(s) at a 
unique point and it is tangent to the graph. 

Then at the moment t = i* the following properties of the function p a (t,x), 
entering the solution (p a ,u a ) to the system (23\) . hold: 
• 

p a (U,x*) ~ Bix^Qfois^a' 22 ^, a -> 0, (44) 

where 

B(x M = k J<^\Z , K m = 1 _ (ml)*r(^ 
\u y (s*)| m 2 2™ my/TT X^™ 

~ means that the quotient of the left-hand side by the right-hand side con- 
verges to 1 as a — > 0; 

p a (t*,x) -» f (s (U,x)), a -> 0, forx^x*, (45) 
where the function sq (t, x) has been introduced in the proof of Proposition 

m 

Proof. Proceeding as in the proof of Proposition [T] we can readily obtain the prop- 
erty (|4"5|) . Thus, let us dwell on the property (jUJ. 

Let us analyze the formula IjlSp at the point (i* , a;*). To that end we note that 
since uq(s*)£* = x* — s* and = — — * s . , for s belonging to the e - neighborhood 
U St (e) of the point s* we have 

1 



uq{s)U + s - x = (ito(s*) + tt' (s*) (s - S*) + Uq (s**) (s - s*) m ) t 

1 



= -- 4 m) («**) (« - s*) m t* - (x - x,), (46) 



nil 
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with s** € U s , (e), e > 0. 
Then from (fT8|) we have 



1 



\J2lxt*0~ 



fo{s)e~ 



(u (.)n.-..) ; 
/o(s)e ^ ds 



( 2,* U-s,) m ) 2 



ds - 



1 



Ms) e 



(-g(')'+'-.) j 



The first integral I\ is equal to 



ds 



h+h. 



where d a (t^,x^) — > B(x*,t*) /o(s*) as a — > 0, where B(a;*,t*) is specified in the 
statement of Theorem [TJ To evaluate Jo- we have used the formula 



/ 



r (-M 

= V f" ; , a = const ^ 0. 



Now let as prove that d^ — > as <r — > 0. 

Let us choose e > so small that for all s € ?7 S , (e) 



C"0 0)*H 



( il (■ 



< ex. 



Then 
\h\< 



1 



^/2nt^a 



l/oOOl 



\Ms)\ e 



(« (.)t + -..)- ! 



+ e 



R\I7». (e) 

The first part in the right-hand side of the inequality due to the boundedness of /o 
is less than const • e, the second part tends to zero as a — > due to the boundness 
of /o- Since e can be chosen arbitrarily small, the statement is proved. □ 

Remark 6. The following asymptotics of K m holds: 



\[2 

K m = m, + 0(lnm) m — > oo. 



Theorem 2. (Amplitude of the S— function) Let the initial data (/□, uq) for 
the system (Q]j, @) be C 1 - smooth and bounded and let the critical instant t* — 

inf ( — ft — y ) be positive and finite. Assume that the initial datum uq is linear on 

the segment Q = (xi,X2), moreover, the second left-hand derivative Uq(x\ — 0) at 
the point x\ and the second right-hand derivative Uq(x2 + 0) at the point X2 do not 
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vanish. Let x* be the unique point such that the line y = s and the graph of the 
initial velocity y = uq(s) have a common linear segment Q = [sijSa]- 

Then at the moment t = t* at the point x — x* the component of the density 
develops a S - singularity of amplitude 

s 2 

A(x*) = ( f{s)ds. (47) 



Proof. We are going to prove that 

' P ^x)^x)dx A { x^), a^O. 



From (|18p we have 



1 f f (^(sjt.+s-x) 2 

p{t*,x)(j>{x)dx = -_ _ / / /o(s)e 4>(x)dsdx 
y/2nt^a J J 

R R 

1 /" /" (u (.)t.+ — )' 

/o(s)e 2 " 2t * <£(x) da; + 



R R\n 
First we analyze I\. 



y / 2irt*o~ 

R n 

1 /" /" (i. (s)t.«+s-x) 2 

/o(s)e - J cj)(x)dsdx = h + I2 



1 („ ( 3 )t»+s-x) 2 

A = ^(x*) / /o(«) / ^j-^ 6 ^ dxds- 



1 f f - ("o(=)*»+ = -^) 2 

/o(s)e 2 - 2 '* ds (0(x) - 0(x*)) dx + 



E4„(e) n 



i»q(»)l.-H-»r 

/o(s)e ^ ds(>(x)-0(x*))(ix = /11+/12+/13, 



t/** (e)RWa;, (e) n 

where [7^, (e) is an e - neighborhood of x*. The integral Jn is equal to A(x*)(/f>(x*), 
with A(x*) specified in the statement of Theorem [21 /12 and /13 tend to zero as 
(j — ► 0, as can be shown in a standard way. 
Further, 

si 

1 f f ( U Q( S )t»+ a -x) 2 

J 2 = ^==^ y y /o(s)e ^ 0(x)dsefe + 

B -00 

oo 

1 f f _ (« (s)*.+s-*> 2 

/o(s)e 2^ 4>{x)dsdx = I 2 \ + J22. 



■5 2 



Let us prove that I2 vanishes as a — > 0. Since /21 and -Z22 can be analyzed similarly, 
we consider only Jgx. 

For s e [7 Sl _o(e), e > from m = 2, we have 

u (s)t* + s-x = iw[) 2) (s„) (s-s*) 2 ^ - (x-x*), 
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with s** e U Sl -o(s). 
Thus, 



I r r _( " i — — (»->i)'-(«-.,)) J 

hi = -==^/ (si) / / e ^ dsdcc- 



-oo 



s l 0)t. 9 

J /■ /■ ( " \ (s-sQ-Qe-x,)) 2 



- DC 



(/o(s) 0(ar) - /o(«i) ^(a:*)) dsdx + 
fo{s) 4>{x) ds dx 



— hll + -^212 + ^213- 

To evaluate hn we use the formula 



(°» 2 -*) 2 , V2 l x 

b 1 ds = e 2b 



where a 7^ 0, b ^ are constants, Bi is the Bessel function of the second type ( [30j ) . 
Thus, 

hu = ff4/ o (»i)0(x,)Af(t„a;,), (48) 

where 

M(t. >jB .) = 1 t )y * £, L = [ e-4 JftB (*\ dx . 

247T2 \u" Si ) 2 J 4 \ z / 

R 

The integral in the expression of L converges, since the integrand is finite at x = 
and decays exponentially at infinity. 

The fact that /212 and hi3 vanish as a — > can be proved routinely. □ 

Remark 7. Under the assumptions of Theorem^ one can similarly show (analo- 
gously to ffllfy) that 

Po-(t*,x*) ~ a™ M(x*,t*,m) f (s*) R a (t*,x*), a — > 0, 

where M(x*, t*, m) is a constant depending only on the properties o/uq and R a {t* , x*) 
fends to #(2; — x*) as cr — ► m 2?'(M). T/iws, t/ie 5 - function, arising from smooth 
initial data without linear segments, has initially a zero amplitude. 



7. The Hugoniot conditions and the spurious pressure 

As follows from the results of Sec. 2, if fpp and upp are smooth, they solve the 
pressureless gas dynamics system. Now we ask the question which system satisfies 
the FP-generalized solution with jumps obtained in Sec. 4. 

The system of conservation laws ([l) , @ implies two Hugoniot conditions that 
should be held on the jumps of the solution [55]. This signifies that the solution 
satisfies the system in the sense of integral identities. If we denote by £ the velocity 
of the jump and by [h{y)\ = h(y + 0) — h(y — 0) the value of the jump, then 
the continuity equation and the momentum conservation give [/] D = [fu] and 
[fu]® = [fu 2 ], respectively. 
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In the case u 2 > the velocity is continuous, therefore the Hugoniot conditions 
hold trivially. 

We should check these conditions for the jumps in the case u 2 < 0. An easy 
computation shows that the first one is satisfied: for the point x\ = u\t we have: 

„ (/i + tiHui + ua) - (2/i + h)ui - (/i + f 2 )u 2 

X) = = Ml, 

-Jl 

and for x 2 = (u\ + u 2 )t: 

„ (2/i + / 2 )m + (/i + / 2 )m 2 - Am 

X> = J — — =U\ + U 2 . 

h + j 2 

However, the second Hugoniot condition does not hold. To understand the reason 
for this let us estimate the integral term in (|2H|) in the case u 2 < as a — > 0: 



(u — u(t, x)) 2 P x (t, x, u) du 



/2vrfcr 



/oO)(«oO) - u{t,x)) 2 ( e - M 'l t X~ x) ) ,ls : 



fo(s){{uo(s)-u F p(t,x))+(u FP (t,x)-u(t,x))y (e ^ 



_ (u (s)t + s-x) J 

/o(s)(mo(«) - UFp(t, x)) 1 (e 



2 (u F p(t, x) — u(t, x)) 



, _(i. (.)H.-i)' 

fa(s)(u (s) - u FP (t,x)) e 



(u FP (t,x) - u(t,x)Y 



his) e 



ds 



h + h + h- 

The integrals I 2 and I3 tend to zero as a — > due to properties of the Riemann 
data since x) — > u F p(t, x) for almost all i6t. Let us estimate I\. 



1 f / (ui t + s-x) 2 

/1 = ;==- / /1 (ui - u F p(t,x)) e ^ ) rf, s 



/ 27Tt(T 



112 I 



27rt(j 



(/i+/2)((«i+U2)-u F p(t,a;)r e 



2 / „ ^ I rfs = 



Thus, 



h 



e 2^ 2 t — e 



/^rfa (2/i + / 2 ) 2 

1 fi{h +h)u% ( _("iy> 2 _((«i+«2)*-*) 2 



/2^rta (2/i + / 2 ) 2 



(2/1 + / 2 ) 



((5(x — (ui + u 2 ) i) — <5(x — u\ t)), a — > 0, 
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in the distributional sense. 

Thus, the integral term corresponds to a spurious pressure p(t, x) between the 
jumps x = [u\ + U2) t and x — u\ t, namely, 

P(t,x) = f i{l + Jff 2 (6(x-(u 1+ u 2 )t) - 6(x- Ul t)), (49) 

see Fig 131 



+ u 2 ) t 



u t 
1 



x 



Figure 3. Spurious pressure, u 2 < 0. 

The Hugoniot condition [fu] D = [ fu 2 +p] is satisfied with this kind of pressure. 
Thus, we get the following theorem. 



Theorem 3. The generalized FP-solution to the Riemann problem \25\l . \26}) with 
constant left-hand and right-hand initial states for the pressureless gas dynamics 
system in the case of a discontinuous velocity (u 2 < 0) solves in fact the gas dy- 
namics system with a pressure defined by fijjty - 

The analogous calculations for the case of rarefaction U2 < show that 


If 2 / 

h = / - — / /1 (ui - UFp(t,x)) e 5^ 



/2nta 



+OO 

(/1 + h) (("l + u 2) - u FP (t, x)) 2 ( e~ " 1 ll*t "J ds 



'2-Kta 



1 ( X.n „ ("1 t-*) 2 , X. 2 „ ll"l+»2)<-») 2 

{ui - -) fie 2 ^t - [{ui +U2) - -) fie 5^ 



V27rtcr V t t 

Here we use the FP-solution (/fp, ufp), obtained in Sec. 4. 

Thus, Ii — > as g — > 0, therefore the integral term vanishes in the case u 2 < 0. 
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8. Sticky particles model vs noninteracting particles 

In our model the particles are allowed to go through the discontinuity as one 
particle does not feel the others. However, in the frame of the sticky particles model 
the particles meeting each other are assumed to stick together on the jump [20] . 
The noninteracting particles model and the sticky particles model are equivalent to 
the same system (JT|) , (|2j) for smooth densities and velocities, however, if the initial 
data have a jump, the solutions behavior differs drastically. 

8.1. Riemann problem with constant states. Nevertheless, we can study the 
solution to the Riemann problem in the case of u 2 < for the sticky particles model, 
too, using the solution obtained above for the non-interacting model. Indeed, the 
jump position Xj(t) is a point between X\{t) = (u± + u 2 )t and x 2 (t) = u\t. The 
mass m(t) accumulates in the jump due to the impenetrability of the discontinuity 
as 

m(t) = ( Xj (t) - (ui + u 2 )t)(h + f 2 )+ 

+[u 1 t-x j {t))f 1 + h =-{(ui+u 2 )(h + f 2 )-u 1 f 1 )t+ 

+ x{t)f 2 + h = -[«/] t + [/] Xj (t) + / 3 , (50) 
where [ ] stands for the jump value, and m(0) = /a, Xj (0) = 0. Further, if we change 
heuristically the overlapped mass between x\ and x 2 to the mass concentrated at 
a point (see FigHJ, then from the condition of equality of momenta in both cases 
we can find the velocity of the point singularity: 

(«i + u 2 )(fi + f»)(xj{t) - (ui + u 2 )t)+ 

+u\ fi (uit - Xj(t)) = 
= -[u 2 .f]t + [uf] Xj {t) = m(t)xj(t). 
Thus, to find the position of the point singularity we get the equation 

(\f\Xj(t) - [uf]t + h)xj{t) = [uf]xj(t) - [u 2 f]t, 
subject to the initial data Xj(0) — 0. The respective solution is 

= r)j ([«/] * - h + \[7i ~2[uf]f 3 t + ([ufY-[f][u*f])t^j , if [/] £ 0, 

(51) 

and 

*>® = *Hf'-hY if|/1 ^ a (52) 

In particular, from the latter formula in the case fs = we get the known expression 
for the velocity of the jump 26J: 

• (A - 2U1 + U2 _ "jg/j + U right 

XjW 2 2 

It can be checked that Xx(t) < Xj(t) < x 2 {t)- The condition expressed by these 
inequalities is equivalent to the Lax stability condition m < ij(t) < u± + u 2 . 

The formulas describing the amplitude of the delta- function in the density com- 
ponent and the singularity position obtained earlier in [19) .[16]. [12] give the same 
result. Moreover, our method allows to find the jump position in a unique way (in 
contrast to the method used in [12]). 
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FIGURE 4. Changing the nonintcracting particles model to the 
sticky particles model. 



It is worth mentioning that the spurious pressure Q490 does not arise in the sticky 
particles model. 

Remark 8. The case of the singular Riemann problem with fa ^ is in fact more 
different from the case of the regular Riemann problem than it seems at first glance. 
Let us begin with the constant initial density ([f] = 0), when the trajectory satisfies 
Q5ty). It is natural to set the amplitude of the initial S - function fa greater or equal 
than zero. Then the trajectory x(t) is continuous since the denominator in i5ty) 
does not vanish. However, if we take fa < 0, then the trajectory goes to infinity at 
the finite moment where m(t) vanishes, and then the trajectory jumps to infinity of 
the other sign. 

Further, if [/] 7^ 0, then we have to use the formula h51\) . It can be checked 
that it is possible to find values u\, 112, fx, fi (for example, u\ = — 1, U2 = —2, f\ = 
2, fi = —1.8) such that the expression under the square root vanishes within a finite 
time t* . However, as follows from i5U\) . \51\) . 

m(t) = yVf ~2[uf]f 3 t + ([ufY-[f]Wf])t\ (53) 

therefore within the same time t* the amplitude of the 6 - function becomes zero. 
Thus, at the moment t* we have to set a new Riemann problem with the jump at 
the point xit*). 

8.2. Riemann problem with non-constant states. Now we extend formulas 
(ISTj) and (|52p to the case of the Riemann problem with non-constant left and right 
states: 

fo(x) = h{x) + 0(x)f 2 (x) + f 3 6(x), (54) 
u (x) = ui(x) + 9(x)u 2 (x), (55) 

where u\{x), 1*2(2;), fi(x), f2{x) are smooth functions, ^3 is a real constant. We 
restrict ourselves to a situation that is quite similar to the case of constant states. 
Namely, we assume that for every i£l and t > the straight line y = ^f^- has at 
most two common points with the graph of the function y — uq(s), moreover, let 
= lim u 2 (x) < and assume that the intersection points s_(t, x) and s+(t,x) 

lie on either side of the origin x — (see Fig|5]). 
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Figure 5. The overlapping domain in the Riemann problem 



So, every point x at the moment t > lies in the overlapping domain 

D = [x-(i),x + (t)], X—{£) = it+t, x+(t) = u-t, 

if the line y — and the graph of y = uq(s) have two points of intersection, 
s_(t, x) < and s+(t,x) > 0. We set u_ = lim u\(x) and u+ = m_ + vf^. 

Further, we assume that for any fixed t the number K of points Xk(t) € D such 
that the straight line and the graph of the initial velocity have a common linear 
segment is finite. Let us again denote by Xj (t) the position of the singularity that 
should change the overlapping domain D in the sticky particles model. Then the 
conservation of mass gives 

Xj{t) U—t 

m(t) = J f(s+(t,x))dx + J f(s_(t,x))dx + J2A k (t,x k (t)) + f 3 , (56) 

u+t Xj (t) k=1 

where Ak(t,Xk(t)) is the amplitude of the 5 - function formed at the point Xk(t), 
calculated by using the formula (|T7]) . 

Further, from the conservation of momentum we have 

Xj(t) u_t 

m(t)xj(t) = J f(s+(t,x))u(s+(t,x))dx + / f(s-(t,x))u(s-(t,x))dx+ (57) 

u+t Xj (t) 
K , 

+ Ak(t,x k (t))—x k (t), 
k=l 
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moreover, the velocity —Xk(t) can be found by the formula 



dt 



— X k {t) = Uo(s*(*,Xfc (*))), 



where s*(t, x~k(t)) is the coordinate of the point where the graph of y = uq(s) and 
the line y = have a common linear segment. If the graph of y = uq(s) does 
not contain any linear segments, then the respective parts in formulas (|56[) and 
(|57|) vanish. Thus, it is sufficient to substitute ([56]) into (|57|) to get the integro- 
differential equation that governs the singularity position. This equation should be 
considered together with the initial condition Xj (0) = 0. 

8.3. Evolution of the singularity formed from smooth data. As we have 
seen in Sec|6j if at a point x and a moment of time t* starting from smooth initial 
data loses its smoothness, then there arises a gradient catastrophe in the velocity 
component (the derivative becomes unbounded), whereas in the density component 
there arises a S - singularity. In the framework of the pressureless gas dynamics 
for t > the S - singularities encompass the overlapping domain D and in the 
overlapping domain the spurious pressure (given by the integral term discussed in 
Sec|7]) appears. In the sticky particles model we have to collapse the overlapping 
domain to the one point, where the whole mass of D accumulates. The position 
of this new singularity should be found based on the conservation of mass and 
momentum. 

For the sake of simplicity we assume that every straight line y = intersects 
the graph of the smooth initial velocity y = Uq(s) at most three times. For every 
fixed x initially the intersection point is unique (0 < t < t*), then at the moment 
t = t* the straight line becomes tangent to the graph of the initial velocity in a 
certain point, and for t > i* we have three intersection points. Our aim is to find 
the position of the singularity and the amplitude of the <5 - function in the density 
component. 

We denote as before by X-(t) and x+(t) the endpoints of the domain D, Xj(t) 
- the position of the new singularity, A(x-(t)) and A(x+(t)) - the amplitude of 6- 
functions produced at the point where the graphs of y — uq(s) and y — x± ^^ s 
have a common linear segment (for fixed t > t*). Further, let s_(t, x), s (t,x), 
s + (t,x) be the subsequent point of intersection (s_ < s < s +-) Let m(t) be the 
amplitude of the S - function in the density component. 

Then, the conservation of mass gives 



(f(s+(t,x)) + f{s (t,x))dx + / (/(s_(t,a;)) + f(s {t,x))dx, 



m(t) = A(x-(t)) + A{x+(t)) + 



(58) 





where A(x±(t)) ca be found by formula (|47|) . 

From the conservation of momentum analogously to (|57|) we have 



m(t)±j{t) = J (f(s + {t,x))u(s + (t,x))+f(s Q (t,x))u(s Q (t,x)))dx+ (59) 
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x + (t) 

J (f(s-(t,x))u(s-(t,x)) + f{s (t,x))u(so(t,x))) dx + 

A(x_(t))u (s_{t,x_(t))) + A(x + (t))u (s + (t,x+(t))). 

Thus, equations ([55]) and the initial conditions m(i*) = A(x*(t*)), Xj(t*) = 
x* define the position and the amplitude of the singularity of the S - function in the 
component of density. Here x* is a point such that y = x *~ s and y = uq{s) have a 
common point s* or a common linear segment [s_, s+] such that the derivatives of 
both functions are equal at s* or on [s_, s+] , is defined in the statement of 

Theorem [5J 

Remark 9. If we want to consider the global evolution of the solution to the Burgers 
equation itself we should set fo = const. The continuity equation plays here an 
auxiliary role. We are not interested in the properties of the density f, which is 
constant everywhere except for domains of vacuum and except for the points giving 
rise to the S - singularity. 

9. Extension to more general scalar conservation law 

Let us consider the following equation 

d t v + (G(v),W)v = 0, (60) 

subject to initial data v(x,0) = vq(x), where v(x,t) — (vi,...,v n ) is a vector- 
function R' i+1 — > K™, G(v) is a non-degenerate differential mapping from K n to 
R n , such that its Jacobian satisfies det dG g} v ^ ^ 0, i,j — 1, ...,n. 
We can multiply |60|) by VGi(v), i = 1, ...,n, to get 

d t G(v) + (G(v),V)G(v)=0. (61) 

Thus, we can introduce a new vectorial variable u = G(v) to reduce the Cauchy 
problem for (fBTj) to ^ with uo(x) = G(v {x)). The stochastic perturbation for (|6"Tj) 
is © with U replaced by G(V). 

Therefore we find the representation of the solution to the stochastically per- 
turbed along the characteristics of equation ([BTjl using the formula (fT5|) with G(vo(x)) 
instead of uq(x). 

Thus, we can apply the results obtained in the previous sections to the inves- 
tigation of the Riemann problem and the arising of singularities for the following 
analogue of the pressureless gas dynamics system: 

gt + div x (gG(v)) = 0, (gG(v)) t + V x {g v <g> G{v)) = 0, (62) 

where g(t, x) is a scalar function that can be interpreted as a density. 

Thus, just as we relate with the non-viscous Burgers equation the system of 
pressureless gas dynamics, so also with the equation (|60[) one can relate the system 

To obtain the solution to the Cauchy problem for (|6"0|) itself we have to consider 
the solution to the Cauchy problem for (|6"2"]) with the data (g , G(vo)), go = const 
and then perform the inverse transform v(t,x) = G _1 (u(t, x)). 
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10. Conclusion 

1. Let us notice that the solution of the Riemann problem for the pressureless 
gas dynamics system obtained in this paper for the ID case satisfies the entropy 
condition 

u(t, x 2 ) - u{t,xi) < 1 

X2 — X\ ~ t ' 

for any sufficiently small x\ and x^ (e.g.[19j) and the balance relations on the jump 
that the definition of solution in the sense of integral identity implies (|16j) are 
satisfied as well. It is known that these conditions do not guarantee uniqueness 
([19j. [12]), However, our solution is unique both in the case of rarefaction and 
compression. In the case of rarefaction it is automatically self-similar (we recall 
that the assumption of the self-similarity implies uniqueness 29 ). In the case 
of contraction the problem of uniqueness was open for the solution to the singular 
Riemann problem, where a non-zero mass is concentrated on the jump at the initial 
time. As was noticed in |12j . for the uniqueness one must prescribe the derivative 
of the amplitude of the 5 - function. In our framework the solution is unique and 
the value of the derivative of the amplitude of the 5 - function follows from the 
expression for the amplitude itself. 

2. In [17] an analog of the system (I19| . ([20]) (without the integral term) in the 
ID case was obtained. Namely, it was proved that for smooth initial data a local 
in time strong solution (p(i, x), u(t, x)) to this system can be constructed by means 
of a nonlinear diffusion process 

t 

X t = X a + J E[(u a (X ))\X s ]ds + aW t , C(X ) = p(0, x), 
o 

so that p(t, x) is the probability density of the diffusion process X t and u(t, x) — 
E[(uo(X ))\X t = x] (with E[., \ ..] standing for conditional expectation and £(X ) 
for the probability density of X ). In fact, this result relates to our Proposition 2, 
since we have shown that the integral term arises only for discontinuous data. 

Further, in [TB] the system (TH?)) , ([2"0"|) (without the integral term) was considered 
in any dimension. In this work there was constructed a global weak solution using 
discrete approximations, and the interaction of particles is given by a sticky particles 
dynamics. 

3. There exist formalisms to represent solutions of parabolic PDE's as the ex- 
pected value of functional of stochastic processes (see e.g. [22], [23], [3], [5] and 
references therein). In particular, in |11] one can find a recent result concerning the 
stochastic formulation of the viscous Burgers equation. An alternative approach to 
the stochastic formulation for a much more wider class of parabolic equations and 
systems can be found in 0]. 

4. We would also like to mention the paper [9], where a numerical method of 
particles for the solution of the pressureless gas dynamics in ID an 2D case has 
been developed. The method is mostly inspired by [20] and in fact in this paper 
the problem of transition from the non-interacting particles to the sticky particle 
model was solved numerically. 
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